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ABSTRACT 

We consider the properties of the warped accretion disc in NGC 4258 which is de- 
hneated by maser emission. We use our analytical models to consider whether the 
disc could be warped by Lense-Thirring precession. We show that such models fit the 
shape of the disc well and we determine the goodness of fit for various combinations 
of the warp radius and the disc and black hole configurations. Though the fits are 
compelling evidence, we note that such a model has implications for the formation 
and longevity of the disc which might be problematic for the current understanding 
of Seyfert Galaxies. 
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1 INTRODUCTION 

NGC 4258 is a bright barred-spir al galaxy at a distance 
of 7.2 Mpc (|Herrnstein et al."l999'). It is in Seyfert 's first 
catalogue of active galaxies (Se yfert 1943) . It has anoma - 
lous spiral arms of Ha emission (|Courtes fc Cruvellieij[l96ll ) 
in the inner regions which are symmetric with respect to 
the nucleus. This emission is probably from shocks formed 
where the matter ejected from the nucleu s meets t he in- 
terstellar medium (|van der Kruit. Port fc Mathewsonll 19721 : 
iFord et al.|[l986l ). 

Masers are the microwave equivalent of lasers. Water 
vapour maser emission at a wavelength of 1.35 cm provides 
information about the accretion discs around highly com- 
pact objects in the centres of active galaxies. Masers are 
bright w here there is no gradient in the bulk line-of-sigh t 
velocity (|Grinin fc Grigor'evI [l983l : I Watson fc Wvldll200d ). 
In a nearly edge-on Keplerian accretion disc the strongest 
features occur along the line to the disc centre where the disc 
is moving perpendicular to the line of sight, these are the 
systemic masers. We also see masers where the plane of the 
sky intersects the disc, these are the high-velocity masers. 

Masers were first dete c ted in the galaxy NGC 4 258 b y 
IClaussen. HeiUgman fc Lol l|l984h and lHenkel et all l|l984h . 
The structure and dynamics of these discs can be accurately 
probed with very long baseline interferometry because the 
emission lines are strong and intrinsically narrow. NGC 4258 
has a set of masers that might trace out a warped accre- 
tion disc which i s seen nearly edge on l|Mivoshi et all 1 19951 : 
iHerrnstein et al] |2005| ). Radio interferometry has shown 
that their rotation is Keplerian (^akai, Inouc & Mivoshl 
1 19931 ) and that the disc extends from 2.8 mas to 8.2 mas 
from the nucleus. 

At its distance of 7.2 Mpc, in the plane of the sky we 
have the scaling Imas = 1.077 x 10^^ cm = 0.0349 pc. The 



Keplerian rotation of the high velocity masers implies an 
enclosed mass of 3.78 X 10^ Ma w ithin 2.8 mas of the nu- 
cleus (Herrnst ein. Greenhill fc Mo ran 1996). This is consis- 
tent with a central supermassive black hole. 

On a scale of milliarcseconds the central core breaks up 
into a jet oriented along the axis of the water maser disc 
IHerrnstein et al.lll997l . Il998h . The northern jet has a flux 
of 3mJy and its mean location is about 0.4 mas north (in 
the plane of the sky) of the implied position of the black hole 
which is the centre of the rotation of the masers. However 
the jet position varies in time. The southern jet has not var- 
ied in flux or position and is located about 1.0 mas south of 
the black hole position. It has a flux of 0.5 mjy. The free-free 
absorption in the masi ng disc causes the differ ence in bright- 
ness between the jets (|Herrnstein et af]|l997h . The position 
angles of the jets are poorly determined. We estimate the 
northern jet is at a posi tion angle of 5° ± 15° measured from 
IHerrnstein et al] l|l998l ). The jets do not appear to be very 
well aligned with each other and the southern jet is at a po- 
sition angle of about 173°. O n larger scales th e jet appears 
to be orientated north-south (|Cecil et al.ll2000l ). 

The systemic masers are about 0.57 mas below the 
disc centre as indicat ed by the high velocity masers 
IHerrnstein et all |2005| ) and at a disc radius of 3.9 mas. 
Hence, the masers at the inner edge of the disc suggest that 
it is inclined at an angle of ^in = sin^^ (0.57/3.9) ~ 8.4° to 
the plane of the sky. There appears to be no structure to the 
vertical positio ns of the masers in t he disc and so the disc is 
probably thin iMivoshi et al.lll995l ). 

Masers operate when the temperature of the surface 
layers in the disc is greater than 300 K but less than 1000 K 
|Moran et al.lll995l ). The outer edge of the disc is where the 
molecular to atomic transition occurs, just outside the outer 
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maser. The reason for the apparent inner edge of the maser 
disc is not yet known. 

The high-velocity masers and are not colinear with 
the systemic masers and have negligible accelerations 
l|Greenhill et al.lll995l ). The radial dependence of declina- 
tion with respect to the systemic velocity of NGC 4258 
is anti-symmetric in the the red and blue-shifted masers 
ijMivoshi et al.lll995l ). It is suggested that the rotation axis 
of t he disc varies with radius by an angle of up to 0.2 radi- 
ans (iHerrnstein. Greenhill fc Moranlll996l l and so the disc is 
warped. Herrnstein et al. ( 20051 ) found that the maser spots 
show a deviation from Keplerian rotation of about 9kms~^. 
They modelled this with a warped Keplerian accretion disc 
with a radial gradient in its inclination of 0.034 mas~^. 

There are several sugge sted explanations for the ori- 
gin of the warp in the disc. ICaproni et al.l ((2006') consid- 
ered some mechanisms for the warping and pr ocession of 
galactic accretion discs. IPapaloizou. Terguem fc LinI 11991) 
showed that it could be produced by a binary companion 
orbiting outside the maser disc. Such a companion would 
need a mass comparable to that of the disc but there is no 
observational evidence for it. A second suggestion is that ra- 
diation pressure from the central black hole pr oduces torque s 
on a slightly warped disc and the warp grows (|Pringlelll99^ . 
However the masing disc is stable a gainst this radiation in- 
stability if ai ^ 0.2 (|Capro ni et al. '2007*1, wh ere ai is the 
viscosity parameter (jShakura fc S unvacv 1973 ) . 

Alternatively, in the absence of other torques, 
ICaproni et al.l |2003) concluded that the warping in the disc 
of NGC 4258 is due to the Bardeen-Petterson effect. If we 
have an accretion disc around a misaligned spinning black 
hole, Lense-Thirring precession drives a warp in the disc. 
The inner parts of the disc are aligned with the black hole 
jBardeen fc P etterson"l975''). 

Caproni e t al. ( 2007) found that the warp radius in the 



disc is comparable to or smaller than the radius of the inner 
masers. They find the timescale of alignment of the sys- 
tem to be a f ew billion years. They based their work on 
the results of IScheuer fc Feileil l| 19961 ) who assumed that 
the surface density is constant. We use the more realis- 
tic power law surface density and viscosities described by 
iMartin. Pringle fc ToutI (|2007l ). We fit these analytical disc 
models, warped by the Lense-Thirring effect to the shape of 
the observed maser distribution. 



2 STEADY STATE DISC MODEL 

In this Section we consider some properties of accretion 
discs. We assume that we have a steady state disc where 
= const and that the surface density is a power law 



(1) 



jShakura fc Sunvaev|[l973l ). To be in steady state, the vis- 
cosity must obey 



E_ 

RoJ ' 



(2) 



where (3, Eq and uio are constants and Ro is some fixed 
radius. 



According to IShakura fc SunvaevI (|l973l ) the thickness 
of the disc in the gas pressure dominated region has 



cx r-^(l - r"2)5 (3) 

where r = R/SRg and Rg = 2GMbh/c'^. This is well ap- 
proximated by H/R = const for 7? > Tig. For NGC 4258 
Tig = 1.05 X 10""^ mas which is much smaller that the in- 
ner edge of the disc at 2.8 mas and so we assume that 
H/R = const in our work. 

2.1 Viscosities 

There are two viscosities: vi corresponds to the azimuthal 
shear (the viscosity normally associated with accretion 
discs) and i/2 corresponds to the vertical shear in the disc 
which smoothes out the twist. The second viscosity acts 
when the disc is non-planar. We assume that the second 
viscosity obeys the same power law as i^i so 



1^2 = 1^20 



(4) 



and 1^20 is a constant. 

We use the a-prescription to describe the viscosities in 
the disc so that 



— aiCsH = aiH fl 



(5) 



where TT is the scale height in the disc and Q is the angular 
velocity and similarly 



1^2 — a2CsH — a2H Q, 



where 



ai 



am — and 02 



■ 020 



(6) 



(7) 



with Qio, a20 and x constant (see below) and Tio is some 
fixed radius which we define in Section 13.11 We gen erally 
take aio = 0.2 and Q20 = 2 (|Lodato fc PringldlioOTl) . The 
angular velocity in the disc is Keplerian so that 



fGMmX 



\ R^ J ' 
Because Cs = HQ the azimuthal shear viscosity is 

'H' " 
R 



(8) 



■ ai 



R^n 



= «,o(GMBHTio)M'f) (I 



(9) 



This can be written in the form of equation ([2} with x — 
P — \ which is a constant and so 

vw = aw{GMwRo)^ ■ (10) 

We can find a similar equation for the vertical shear viscos- 
ity, V2- 



2.2 Surface Density 

We let the surface density be the power law given by equa- 
tion ([T]) but we need to find Eq. We have the equation de- 
scribing the hydrostatic equilibrium in the disc 

d7 = -^5- (11) 
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where P is the pressure and is the z component of the 
gravitational force, p is the density and the equation of state 
is 



(12) 



in the case of a vertically isothermal model with the, 
constant in z, sound speed. We find the z-component of the 
gravitational force to be 

GMbhz GMbhz 



for z <^ R. We can integrate equation pip over z using 
equations p2|l and (|13l) to find the pressure 



and density 



P = Po exp 



po exp 



GA/bh 



2 C2J?3 



z^ GMbh 



cms 



(14) 



(15) 



where Po and po are the pressure and density at the mid- 
plane z = 0. We can now find the surface density 



r 



because 



H 



pdz = Vii^Hpo 



, GA/bh, 

For a steady state accretion disc we have 
M = 37ri/iE, 



(16) 



(17) 



(18) 



with j/i the usual viscosity. We eliminate E from equa- 
tions Uni) and (O to find 



M 



GA/bh M 



l|Neufeld fc MalonevlllOOSI ) and, using equation (|16p , we find 



E = 



M 



R 



If H/R — const then this is equivalent to 



where _Ro is some fixed radius and 



Af 



R 



Sttoio \HJ (GA/bh Po) 2 



(20) 



(21) 



(22) 



We have viT, — const which is necessary for a steady state 
disc. 



2.3 Inclination of the Disc 

Following iMartin. Pringle fc ToutI l|2007l ) we consider the 
disc to be made up of annuli of width dR and mass 27rEPdP 
at radius R from the central object of mass M-qb with sur- 
face density Ti{R,t) at time t and with angular momentum 
L = (GA/bhP)^''^EZ = LI. The unit vector describing the 
direction of the angular momentum of a disc annulus is given 
by / = {lx,ly, h) with |/| = 1. 



We use equation (2.8) of lPringld (|l992l ) setting dL/dt = 
and adding a term to describe the Lense-Thirring preces- 
sion (the last one) to give 



0=1 A 

RdR 



3P 9 , ^, 3 



T 1 r,r dl 



+ 



LJr, X L 



P3 



The Lense-Thirring precession is given by 

2GJ 



(23) 



(24) 



iKumar fc Pringld [l985l ) . where the angular momentum of 
the black hole J = Jj with j = {jx,jy,jz) and Ij'j = 1 can 
be expressed in terms of the dimensionless spin parameter 
a such that 



J = acA/BH (^f^) 



(25) 



We use the frame of the black hole with the black hole 
at the origin and with i ts spin along the z axis s o that J = 
(0, 0, J). In this frame, iMartin. Pringle fc ToutI (|2007l ') find 
the steady state solution of equation (|23|) for W = Ix + Hy, 
with i = V—T, to be 



W{R) = 



2Wo. 



(_i)4(l+;3) 



„ ( V2 R 
X K 1 -(1 — i) — 

2(TW I 1 + /3^ ' \Ro 



1/4 



(26) 



where Waa is the constant value of as P — > cxa. It rep- 
resents the direction of the angular momentum of the outer 
disc. This solution has the outer disc angular momentum 
vector in the x-z plane. Here Ki, is the modified Bessel func- 
tion of order v and F is the gamma function. In deriving this 
solution we assume that we can neglect the non-linear term 
l.d^l/dR^ — ~\dl/dR\^ . We need the warping to be gradual 
enough for this to be true. 

We let the outer disc have the fixed angular momentum 
direction vector {IxocOJzoo)- As P — > oo, W ^ Woo ~ 
Ixoo ~ sin rj where rj is the inclination of the outer disc to 
the black hole. The inclination of the disc to the black hole 
at radius P (the z-axis) is 



e{R) = cos"'(/,) = cos"'(/r~fWp). 



(27) 



We let the angle <j> be the usual azimuthal angle in spherical 
polar coordinates with ^ <^ < 27r. It is the angle in the x-y 
plane from the projected direction of the outer disc angular 
momentum which in our model, given by equation (|26p . is 
along the a;-axis. 



2.4 Disc Section 

The high velocity masers trace the disc in NGC 4258 in 
the plane of the sky so we need to find the cross section of 
our disc model in this plane. Our solution is in coordinates 
{x,y,z) in the frame of the black hole and the outer disc lies 
in the x-z plane. We initially take this to be aligned with 
the coordinates {xhh,ybh,Zhh) which are the coordinates that 
we take the disc section in the Xbh-Zbh plane. 

We rotate the coordinates (x,y,z) about the Zbh axis by 
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an angle 
given by 



where . 



The shape of the disc in the a:;bh-2bh plane is 



R 

Zh\i[R) = —-5- tan^ cos( 



0n 



(28) 



= arg(VK) and we have the negative sign so that 
our model has the same orientation as the maser data and 



2;bh(-R) = i5-(sin^( 



- 0max) + COS^ 6lcOS^((?ic " 'i/lmax))^ 

(29) 

in the positive Xhii direction. We are assuming that the disc 
is symmetric so the side of the disc with negative Xhh can 
be found by rotating this shape by 180° about the j/bh axis. 

\i 4>c = Q we see the largest warp possible for a given 
disc because we are taking a disc section in the plane which 
contains both the black hole and the outer disc angular mo- 
mentum. Conversely, if (j>c — 'k/2 then we see the least warp- 
ing in our disc section. 

We next consider the disc section when the black hole 
is not in the plane of the sky but inclined at angle iinc to 
it. We rotate the {x, y, z) coordinates about the Xbh axis 
through an angle —iinc. For small iinc the shape of our disc 
becomes 



R 

ZhtiiR) = - — tan^ 
Ho 



nax ) cos ^inc 



(30) 



and ajbh remains the same. If the plane of the sky contains 
the black hole spin then iinc = 0. We can use this small 
angle approximation because if the angle iinc were big then 
we wouldn't see the masers because they are only seen when 
there is no gradient in the bulk line of sight velocity. This 
only happens when the disc is close to edge on. 

We sketch the plane of the sky coordinates (a;bii, Zbh) 
and the coordinate system for the disc {x, y, z) in Figure [1] 
We show the plane of the sky, P, in which we find the disc 
section. It contains the dashed lines. The normal, n, to the 
plane P is the line of sight. The disc section is in coordinates 
2;bii and ^bii- 

In the case where iinc ~ and 0c — 0, we find the 
inclination of the disc at radius R to the plane of the sky to 
be 



( = n/2 — cos ^{ly)- 



(31) 



The observations of the systemic masers predict that, at 
R = 3.9 mas, the disc has inclination ("in = 8.4°. For our 
models with iinc ~ and 0c = we find the inclination of 
our disc at R — 3.9 mas and compare this to the observed 
value. 



3 PROPERTIES OF NCG 4258 

The inner radius of the masing disc correspo nds to Rjn — 
2.8 ma s and the outer edge to i?out = 8.2 mas. lArgon et al.l 
l|2007l l examine the systemic masers in the direction per- 
pendicular to the axis of the disc. They observe that the 
vertical thickness of the disc has a limit H < 0.005 mas. Be- 
cause the systemic masers are at i? = 3.9 mas we see there 
is an observed upper limit on H/R such that 



R 



< 0.0013. 



We use a canonical value of H/R 
calculations. 



0.001 in the following 




n 



Figure 1. A diagram showing the plane of the sky, P, with nor- 
mal n which corresponds to the line of sight. We find the disc sec- 
tion in the plane P. The Cartesian coordinates used for the gen- 
eral solution to the warped disc equation given by equation II26II 
are (x, y, z). The black hole spin is up the z axis and the outer 
disc angular momentum is in the x—z plane. The coordintates in 
the plane of the sky are Xbh and Zbhi and j/bli is in the direction 
of —n. To find the orientation of the black hole and the disc we 
start with {x,y,z) aligned with (aibii, yblu -^bli)- We then rotate x 
and y about the Zbii axis by an angle —0c. Then we rotate about 
the a'bii axis by an angle -iinc- The dashed lines are in the plane, 
P. 



3.1 Warp Radius 

The warp radius is where the Bardeen-Petterson effect is 
balanced by the viscous evolution of the disc and is given by 

L0„ 



Rwa,rp — 2 (*^^) 
1^2 (-ftwarp ) 

(|Scheuer fc Feilerl Il996l : iMartin. Pringle fc Tout! |2007| ). 



(34) 



tion ^ for the second viscosity we find 



and, rearranging, 

-^wari 



4aG^MgH fRV /jRLtp \ 



C-a20i?larp \HJ \GMbhJ 



ia{GMB 



More specifically, for NGC 4258 we find 



-f^warp — 0.83 



(35) 



(36) 



This warp radius is closer to the black hole than the inner 
masers and we consider this in more detail in Section [3.41 
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3.2 Accretion rate 

The motion of the infalhng gas in NGC 4258 cannot be 
measured directly so we have to estimate the mass accretion 
rate, M, indirectly. We use the bolometric luminosity of the 
system and assume that 



Lboi = eMc 



(37) 



where e is the accretion efficiency of a Kerr black hole, e ~ 
0.1. The bolometric luminosity of NGC 4258 is about 4 x 
10*'^ ergs~^. This corresponds to a mass accretion rate of 

M = 7 X 10"^ Mq yr"^ (38) 

if the disc is in a steady state. 

3.3 Surface Density and Mass of the Disc 

We explore the consequences of the assumption that the 
disc is in steady state. When we vary /3 we want properties 
of the disc at i? = -Rwarp to be be fixed and so we choose 
Ro = -Rwarp . We find the surface density using equation (|2ip 
to be 



E = Eo 



R 



R^ 



(39) 



where 



M f R 



37rai \H) {GMmR^.r^)^ 



= 1 



V 0.2 / V 2 / VO-OOl / 



M 



gem 



Jx lO-5M0yr-i/ ■ ^^^^ 

ICaproni et al] (|2007l ) used the same model for the disc with 
Qi = const but varied H/ R with a power law in R. 

Using the surface density of the disc we can find the 
total mass of the observable disc to be 

Md = 27r y^"""" ^RdR = ^i?earp(i?o:/ - I^-^) (41) 

if /? / 2, or, if /3 = 2 then 



Md = 2^Eoi?Lrp log 
For NGC 4258, with /3 = 2, we find 



(42) 



Md =2.97 



H/R 
0.001 



M 



7 X lO-5M0yr-i 



(43) 



This is much smaller than the mass of the black hole but 
in the next section we consider if this is, or could be, large 
enough for our disc to be self-gravitating. 



3.4 Self Gravity of the Disc 

Having calculated the mass of the observed maser disc we 
check if the disc is stable against the effects of self gravity. 
If it were unstable against self gravity then it would collapse 



into bound fragments locally. If the mass of the disc satisfies 
the inequality 



Md < Mbh ( ^ 



(44) 



then the disc is stable (jToomrd [l963 : iBinnev fc Tremaind 
198'J). For H/R = 0.001 this corresponds to a critical disc 
mass of 3.78 x IO^Mq. For /3 = 2, using equation (|43l), we 
see that we need H/R > 0.0006 for a disc which is stable 
against self gravity. Now we have constrained H/R to the 
approximate limits of 



0.0006 < H/R < 0.0013 



(45) 



if our model is correct. Using equation (|36p we find that for 
13 = 2 our warp radius is limited to 



0.58 < 



1 mas 



< 1.65. 



(46) 



This is a relatively small range. We note that this range for 
flwarp is comp l etely i nside the inner observed maser points. 
ICaproni et al.l l|2007h found that the warp radius was ap- 
proximately equal to or smaller than the inner maser but 
we find that it is at most around half way between the black 
hole and the inner maser. 

We note that this mass just takes into account the mass 
of the disc of the observed masers. For our steady state 
solution, the surface density increases as we move towards 
the black hole and so, if we calculate the mass of the disc 
including inner regions, then our disc soon becomes self- 
gravitating. If the inner parts of the disc were not in steady 
state or if H/R were much larger there then we could over 
come this problem. 



3.5 Gravitational Torques 

In our model, different radii of the disc communicate 
through the two viscosities. Discs of non-negligible mass 
could also communicate through gravitational torques. We 
consider he r e how these torques would affect equation (|23|) . 
iNavakshinI l|2005h considered the gravitational torque be- 
tween two rings at radii R and R\ of masses M and M\ 
respectively which are inclined at an angle 7 to each other. 
The torque exerted by the second ring on the first is 



47J-2 



ri X r 

\r — ri\^ 



d(j> d(j)i 



(47) 



where r and ri are the position vectors around the rings 
and the integration goes over and (j}i which are azimuthal 
angles in the frame of the respective rings. He found the 
maximum precession frequency to be when R/Ri = \/3/7 
and the maximum torque to be 



Tgrav = 0.085 Am^R^ sin 7 cos 7 ^'^^ 



0.085 Mf7^i?^ sin 7 cos 7 



Mbh 
27ri?iEi(5i?i 



(48) 



where Mi = 27rEi J?i(5_Ri is the mass of an inclined annulus 
on which we consider the torque and Ei = E(J?i). In this 
formula M is the mass of the annulus at R which provides 
the greatest torque. We replace M by the mass of the disc, 
Md to obtain an upper limit on the total torque on the 
annulus at R\. 
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The magnitude of the viscous torque is (equation (f23l 

dl 



T'. 

VIS 



1 1 
1 



dR 



2RI 



,2 V2iRi)Li sin 7 



because 



Ri 



dl 



OR 



■ smj. 



(49) 



(50) 



We multiply by the area of the ring to find the viscous torque 
acting on a ring of width 5Ri at radius Ri to be 



2nRiSRi^,y2{Ri)Lism-y. (51) 



We now compare the magnitude of the maximum gravita- 
tional torque to the viscous torque 

T 0. 085 A/d -R^ i^i sin 7 cos 7 ^ 1 

Sgrav — rr — f D \ T ■ "> ^ 

ivisc Tvu2(Ri)Li am-y 

and subsituting we find 



f Tgra,v 

?grav „ 

J vise 



0.085 MM^R^Ri sin 7 cos 72ig|i 



7r!/2(-Ri)(GA/BH-Ri)2 El sin 7 
0.17 M^GRj cos 7 



U2{Ri)R{GMBnRiy- 

0.17 M^Rf cos 7 
i^2oR{G-^MbhRi)^ \R- 



Ri 



(53) 



With (3 = 2 and cos 7 ~ 1 at i? = Rin and Ri = 
^/T/SRin we find that ^g^av = 3 x lO"'*(Md/M0) and at 
= 7?„^t and R = 'jtj7Ri then ^g^av = lO"''(Afd/M0). 
Using the disc mass derived in Section 13.31 of A/d = 3 x 
10^ M0 we find at i? = _Rin, ^grav = 0.9 and at R\ = Rout, 
^grav = 0.3. This disc mass is the critical one at which the 
gravity torques become as important as the viscous torques. 
For smaller disc masses self gravity torques would be unim- 
portant. In any case the gravitational torque is not as large 
as the viscous torque even at its maximum so we are justi- 
fied in neglecting it. However, it ought to be included in any 
numerical models especially for larger disc masses. 



4 COMPARING THE MODEL TO 
OBSERVATIONS 

In this section we compare the maser data to our analytical 
disc models and find parameters which give the best fits. 



4.1 Maser Data 

We use th e maser distr i bution data tabulated in the online 
version of lArgon et al.l (|2007l ). This lists positions of each 
maser in the plane of the sky, Zm and Xm. For the first set of 
data, BM056C, we correct the north-south position of the 
masers (Alice Argon, private communication). They corre- 
lated the first three epochs at a less accurate position than 
later epochs and the data in the online table for BM056C 
refiect this. From the relativistic velocity, Urd, in column 3 
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Figure 2. The maser distribution in space. The top plot shows 
all the data points and the bottom plot shows the binned data 
with error bars. 



of Table 5 in the online table of lArgon et al.l l|2007h . we find 
the frequency 

/froq = 22.23508 ( ^ ~ •"rci/e ^ ^ ^^^^ 



where c = 2.998 x 10^" cms^^ . Then we apply the correction 

1^ mas (55) 



Zm — 66.0 



22.19856 



and remains the same. 

In order to get an estimate of the measurement error 
we bin the maser data points. We put the data into 11 bins 
and find the average Xi and Zi in each where i = 1,2. ..11. 
We then find the standard deviations of each average point, 
aix and aiz . We list these in Table [T] In Figure [5] we plot 
all of the maser data in the top plot and below we plot the 
binned data points with error bars. 

The symmetry of the velocity distribution of the masers 
tells us that the black hole is on the line Xn^ = and we 
find the position of the black hole, Zadd, from the symmetry 
of the shape of the maser distribution. We use the inner 
three points on the right in the binned data and find an 
interpolated corresponding Zm value on the left hand side at 
the positive Xn^. We can then find the midpoint of these two 
Zm values from the left and right hand sides. We average 
over the three points calculated and find that Zadd = 0.52 in 
the frame of the masers. 



4.2 Disc Section Transformation 

We have the shape of the disc section as given in the co- 
ordinates Zhh and Xhh. However, we can vary the size, the 
rotation angle and position within the plane to get a good 



find the disc section we want. 

We rotate this disc section by an angle q about the line 
of sight. The black hole spin was originally in the north- 
south direction but with this transformation it is at a po- 
sition angle of —q. If we assume that the observed jets are 
parallel to the spin of the black hole, then the black hole 
spin should be at the same position angle as the jets. 
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Table 1. The binned maser points with the number of points in each bin and the standard deviations. 



We can choose the factor / = i?warp/nias which deter- 
mines where /?warp is relative to the maser distribution. Our 
model becomes 

zt = ±(2:bh cos q + a^bh sin </)/ -I- Zadd (56) 

and 

xt = ±(2;bh cos q - Zbh sin g)/, (57) 

where the sign depends on which side of the disc we are 
modelling. 

We now have a model which can be rotated, stretched 
and moved to fit the masers. We choose the parameters 77, 
/ and /3 for the disc we wish to fit. Then we choose a disc 
section by choosing and iinc- We then rotate the disc 
section to the best fit by varying q. 

4.3 Comparison of our Model to the Maser Data 

We can compare our analytical model to the binned points 
by finding 

i — 1 * 

where = A'^ — m, A'' = 11 is the number of bins and m is 
the number of parameters to be determined from the data. 
We can change the inclination of the black hole relative to 
the outer disc, r\, and also the power /3. We can rotate the 
disc about the axis of the black hole by changing angle ^c. 
We can also vary the inclination, iinc, of the black hole to 
the line of sight. We can vary the warp radius through / . 

For a given set of parameters for 77, /3, /, iinc, ^ic, we 
minimise Xio/10 by varying q. We vary one parameter at 
a time and so m = 1 and v — 10. The critical probability 
values, p, for deviations of 1, 2 and 3 a for Xio/10 given 
in Table [2] The probability of exceeding the critical value of 
Xio/10 is 1 - p. 

In Figure|3]we plot contours of Xio/10 in / and rj space 
for different values of /3, and i. Each point in the plot is 
the best fitting point when we vary q. The plots would be 
symmetric over the line rj = 90° if we were to extend them 
to higher r) because the disc shape is the same for t] and 
■K — Tj but if 77 > 7r/2 the disc is counter aligned. We see that 
we can fit the disc shape for a variety of parameters. 

Generally, we see that we need rj > 30° or so for a 
reasonable fit and the warp radius must be on the scale of 
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0.683 
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0.955 


1.865 
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0.997 
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Table 2. The probability of exceeding the critical value of Xio/^O 
is 1 — p. We tabulate critical values of Xio/10 deviations of 1, 
2 and 3 a. 

the masers. Comparing the top plots which have different 
P only, we see that for the lower value of /? = 1 we can fit 
higher Tiwarp better but we cannot fit such small 77 as for 
/3 = 2. 

The bottom left plot is the same as the top right plot 
except the black hole is inclined at an angle of iinc = 10° to 
the plane of the sky. Comparing this to the top right plot 
which has iinc = we see that the the range of parameters 
we can fit is not changed much for small inclinations of the 
black hole from the plane of the sky. Note that the method 
we have used to generate these fits does not take into account 
whether masers would actually be seen for the inclinations 
chosen. We need the inclination of the black hole from the 
plane of the sky to be small in order to see masers at all. 

In all the plots considered so far, we have chosen the 
disc section of a given disc with the most warping, the one 
where the plane of the sky contains the angular momentum 
of the black hole and outer disc vectors. In the bottom right 
plot we choose a different disc section with (jic = 45°. Here, 
we cannot fit such small J?warp or r). 

In Section 13.11 we found an upper limit on the warp 
radius for the disc to be stable against self gravity that JJwarp 
must be less that 1.65 mas. This rules out the case with 
(f)c = 45° completely. We find that we can only fit discs 
which have |(j!)c| < 25°. In order for our models to explain 
the warping in NGC 4258, the black hole and outer disc 
angular momentum vectors must be close to lying in the 
the plane of the sky. By the bottom left plot, the limit also 
implies that we need 77 > 60°. 

4.4 Best Fitting Models 

We consider the best fitting combination of parameters for 
the case (3 = 2, i = and (f)^ ~ 0. The parameters giving the 
smallest Xio/lO of 0.173 are 77 = 45.8° and i?warp = 5.8mas 
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Figure 3. Contour plots of the values of Xio/W for combinations of i?warp/mas and rj, the angle between the black hole and outer disc 
angular momenta. The solid lines are Xio/10 = 1.154, the dotted lines are Xio/10 = 1.865 and the dashed lines are ~ 2.661 

which correspond to 1, 2 and 3cr deviations. The top left plot has f3 = 1, (f>c = and ijnc = 0. The top right plot has [3 = 2, <j)c = and 
'inc = 0. The bottom left plot has /3 = 2, 0c = and line = 10°. The bottom right plot has (3 = 2, 4>c = 45° and = 0. The vertical 
long dashed lines shown for the cases with (3 = 2 are the upper limit on i?warp derived from the condition that the disc must be stable 
against self gravity. 



and this has q = 10.6°. The position angle of the black 
hole is —q and so if we assume that the jets are parallel 
to the spin of the black hole then the position angle of the 
jets here is —10.6°. This is on the edge of our error for the 
observed position angle of the jets. The inclination of the 
disc a.t R = 3.9 mas to the plane of the sky is ("in = 7.4°. 
In the top of Figure [4] we plot this fit with the maser data 
points. For this warp radius we would need H/R = 0.0002 
and so the disc would be unstable against self gravity. 

The best fitting model which is stable against self grav- 
ity has iiwarp = 1.3 mas, 77 = 85.9°, Xio/10 = 0.20331 and 
q — 39.6° and we plot this in the bottom plot in Figure [l] If 



we take into account our constraint on -Rwarp then we need a 
larger angle q which means that the black hole has position 
angle which is more negative and not within the error on the 
observed position angle of the jets. However, the inclination 
of the disc at i? = 7 mas to the plane of the sky is ("in = 8.6° 
which is closer to the observed value of 8.4°. 

The angles of the jets are very poorly determined and 
vary in time indicating that the angle is affected by the 
surrounding medium. None of the models that fit the data 
well have negative q. Once we constrain the warp radius we 
find we need an even larger q to fit the data. However, the 
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Figure 4. The maser points with the best fitting models with 
13 = 2, line = and ct>c = 0. The top plot has xlo/^^ = 0.173, 
rj = 45.8°, -Rwarp = 5.8 mas and q = 10.6°. The bottom plot is 
the best fitting model which is stable against self gravity and has 
xlo/W = 0.203, r; = 90°, i<;„arp = 1.3 mas and q = 39.8°. 



jet is not very clearly observed (jHerrnstein et al]|l997^ and 
so we cannot rule out our models because of this. 



5 PHYSICAL IMPLICATIONS FOR THE DISC 

We now consider the physical implications of our fit to the 
disc structure in NGC 4258. First we calculate the viscous 
timescale in the disc. For -Rwarp ~ 0.83 mas with H/R — 
0.001 we find the viscosity to be 

ui = 4.3 X 10'** (-^^y cras-\ (59) 

The viscous timescale in the disc for /3 = 2 is independent 
of R and is given by 



= 6.3 X lO'^yr. 



(60) 



and represents the timescale of the movement of mass 
through the disc. The second viscosity timescale is 



— — 6.3 X 10 yr 

Ce20 



(61) 



and this is the timescale associated with the generation of 
the warp in the disc. 

The precession timescale due to the Lense-Thirring ef- 
fect is given by 

(62) 



2GJ 2aG2 M| jj ' 



At -Rwarp we have tur — ti,^ because of the definition of 
^warp. At -Rin, tur = 4.6 X 10* yr which is longer than f^j. 
Hence, for the disc to be in a steady state we require it to 
last for much longer than the viscous timescale. This means 
that the disc must be effectively of infinite mass because it 
is continually supplied with matter. 

Seyfert outbursts are stochastic randomly-orientated, 
short-lived accretion events. If this system is a Seyfert sys- 
tem then the maser disc is just the outer regions of a disc 
generated by an accretion event. Seyfert episodes thought 



to be shorter than 10^ yr l|Sanderslll984h and nuclei evolve 
through at least 100 of these. The disc becomes stable and 
thin in a few orbital periods, a time of around 10^ yr at the 
outer edge of the disc. The warp then just represents the dis- 
tribution of angular momentum in the accreted matter. The 
time to change the warp, around 10* yr, is much longer than 
the length of the Seyfert outburst and the typical time to 
the next outburst. It is hard to reconcile our model without 
upsetting this understanding of Seyfert outbursts. 

However NGC4258 is a rather inactive AGN. The as- 
sumed luminosity is well below the Eddington luminosity. 
We suggest that it could be in a more quiescent phase with 
steady accretion. 



6 CONCLUSIONS 

We can fit the shape of the warped disc in NGC 4258 with 
steady state accretion disc which has a misaligned spinning 
black hole that causes the disc to warp by Lense-Thirring 
precession. We need the angle between the black hole and 
the outer disc to be greater than 60° in order to get a good 
fit. We need the warp radius to be 1 < -Rwarp/mas < 10. 
We need the angular momentum of the outer disc to be at 
an angles of less than 25° to the plane of the sky. The black 
hole spin must be close to the plane of the sky to see masers. 

The best fitting model has the warp radius at 5.8 mas 
and a black hole position angle of —10° but this disc would 
be unstable to self gravity. The best fitting model which 
would not be self gravitating has -Rwarp = 1.3 mas but the 
position angle of the black hole is —40° . The observed direc- 
tion of the jet is 5° ± 15° and so cannot be easily accounted 
for using this model. 

We need the disc to be very long lived in order for it 
to become warped in its lifetime. If the disc was in steady 
state all the way in, it would be unstable against self-gravity. 
The requirements of the model are hard to fit with our un- 
derstanding of Seyfert outbursts which are thought to last 
only 10^ yr. However, it remains that we can fit the observed 
maser distribution with an accretion disc warped by Lense- 
Thirring precession remarkably well. 
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